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Thermodynamic observables of Mni2-acetate calculated for the full spin-Hamiltonian 
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35 years after its synthesis magnetic observables are calculated for the first time for the molecular 
nanomagnet Mni 2 -acetate using a spin-Hamiltonian that contains all spins. Starting from a very 
advanced DFT parameterization [Phys. Rev. B 89, 214422] we evaluate magnetization and specific 
heat for this anisotropic system of 12 manganese ions with a staggering Hilbert space dimension of 
100,000,000 using the Finite-Temperature Lanczos Method. 

PACS numbers: 75.10.Jm,75.50.Xx,75.40.Mg 


I. INTRODUCTION 

Density Functional Theory (DFT) has greatly ad¬ 
vanced over the past years and is nowadays able to 
predict parameters of spin-Hamiltonians with which the 
low-temperature physics of correlated magnetic materi¬ 
als can be described, compare e.g. Refs. 1-15. Along 
this line the complex spin-Hamiltonian of one of the 
most exciting magnetic molecules, Mni 2 -acetate, was re¬ 
cently predicted.^® These calculations consider almost 
all terms that are bilinear in spin operators as there 
are the Heisenberg exchange interaction, the anisotropic 
antisymmetric exchange interaction and the single-ion 
anisotropy tensors. The new calculations outperform 
earlier attempts^^’^® and provide much richer electronic 
insight than parameterizations obtained from fits to 
magnetic observables^® or parameterizations resting on 
knowledge from similar but smaller systems. But despite 
all the success DFT is not capable of evaluating mag¬ 
netic observables which is the reason for the detour via 
spin-Hamiltonians. 

The magnetism of anisotropic molecular spin systems 
is fascinating due to interesting phenomena such as bista¬ 
bility and quantum tunneling of the magnetization.®® 
Bistability in connection with a small tunneling rate leads 
to a magnetic hysteresis of molecular origin in these 
systems. That’s why such molecules are termed Sin¬ 
gle Molecule Magnets (SMM); Mni 2 -acetate is the most 
prominent SMM.®^“®^ But although Mni 2 -acetate con¬ 
tains only four Mn^^^^ ions with s = 3/2 and eight Mn^^^^^ 
ions with s = 2, it constitutes a massive challenge for the¬ 
oretical calculations in terms of spin-Hamiltonians since 
the underlying Hilbert space of dimension 100,000,000 
is orders of magnitude too big for an exact and com¬ 
plete matrix diagonalization.®® But thanks to the fact 
that the zero-field split ground-state multiplet is ener¬ 
getically separated from higher-lying levels a description 
using only the S' = 10 ground-state manifold is suffi¬ 
cient to explain observables at low temperature - this ap¬ 
proach was used in the past. Thermodynamic functions 
which involve higher-lying levels, for instance observables 
at higher temperature, can of course not be evaluated in 
such an approximation. 

Fortunately not only on the side of DFT progress has 


been made in past years, but also in terms of pow¬ 
erful approximations for spin-Hamiltonian calculations. 
For not too big systems with Hilbert spaces with di¬ 
mensions of up to 10^® Krylov space methods such as 
the Finite-Temperature Lanczos Method (FTLM) have 
proven to provide astonishingly accurate approximations 
of magnetic observables.®®“^^ While FTLM has been used 
for Heisenberg spin systems mostly, very recently the 
method was advanced to anisotropic spin systems.^® 

In this article we therefore employ the most recent 
FTLM in order to study thermodynamic functions of 
Mni 2 -acetate starting from parameterizations provided 
by DFT or other methods. We evaluate the magneti¬ 
zation as well as the specific heat both as function of 
temperature and field and compare the various parame¬ 
terizations of the spin-Hamiltonian. 

The article is organized as follows. In Section H the 
employed Hamiltonian as well as basics of the Finite- 
Temperature Lanczos Method are introduced. Sec¬ 
tion HI, IV and V discuss the effective magnetic mo¬ 
ment, the magnetization and the specific heat, respec¬ 
tively. The article closes with summary and outlook. 


II. FTLM FOR ANISOTROPIC SPIN SYSTEMS 

For Mni 2 -acetate, which is a highly anisotropic spin 
system, the complete Hamiltonian of the spin system is 
given by the exchange term, the single-ion anisotropy, 
and the Zeeman term, i. e. 

= X 5 * ■ + X 

i<j i 

+PbB g,sl . 

i 

Jij is a 3 X 3 matrix for each interacting pair of spins 
at sites i and j which contains the isotropic Heisenberg 
exchange parameters together with the anisotropic sym¬ 
metric and antisymmetric terms. In the sign conven¬ 
tion of (I) a positive Heisenberg exchange corresponds 
to an antiferromagnetic interaction and a negative one 
to a ferromagnetic interaction. denotes the single-ion 
anisotropy tensor at site i, which in its eigensystem e/. 
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ef, ef, can be decomposed as 

Di = Aef 0 ef + {e/ 0 e/ - 0 e^] . (2) 

The terms gi could in general be 3 x 3 matrices, too, but 
for the sake of simplicity it is assumed that the gi are 
numbers and moreover that gi = 2 for all ions. This as¬ 
sumption is justified for the and ions in 

Mni 2 -acetate, since the g-factors of both ions are esti¬ 
mated to be very close to 

The Finite-Temperature Lanczos Method (FTLM) ap¬ 
proximates the partition function in two ways:^®’^° 

Z{T,B)^'^f: E e-'’‘-'’K»(>')k>t. (3) 

1/—1 n —1 

The sum over a complete set of vectors is replaced by a 
much smaller sum over R random vectors \v). The ex¬ 
ponential of the Hamiltonian is then approximated by its 
spectral representation in a Krylov space spanned by the 
Nl Lanczos vectors starting from the respective random 
vector \v). \ n{v)) is the n-th eigenvector of H in this 
Krylov space. 

It turns out that a very good accuracy can already be 
achieved for parameters i? « 10 and Nl ~ 100 , especially 
in cases when the low-lying energy spectrum is dense. 


III. EFFECTIVE MAGNETIC MOMENT AS 
FUNCTION OF TEMPERATURE 

Mni 2 -acetate contains four ions with s = 3/2 

and eight ions with s = 2. Following Ref. 16 the 

ions and the exchange pathways are depicted in Fig. 1. 
Mn(iv) 

ions (1-4) are shown as red circles, ions 

(5-12) as blue ones. An S '4 symmetry of the molecule is 
assumed. 



FIG. 1. (Color online) Schematic structure of Mnl2, same 
labeling as in Ref. 16. ions (1-4) are shown as red 

circles, ions (5-12) as blue ones. An 54 symmetry is 

assumed. 


Since the discovery of the pronounced SMM-properties 
of Mni 2 -acetate several groups developed parameteriza- 
tions of the full spin-Hamiltonian. These data sets, of 
which the most prominent ones are given in Table I, con¬ 
tain parameterizations of Heisenberg models and were 
put forward following various scientific reasonings. Ear¬ 
lier attempts assigned values of exchange interactions 
in analogy to smaller compounds with similar chemical 
bridges between the manganese ions. Later investiga¬ 
tions combined for instance high-temperature series ex¬ 
pansion with the evaluation of low-lying excitations seen 
in Inelastic Neutron Scattering (INS) experiments.^® A 
necessary condition that has to be met by all parame¬ 
terizations, is that the ground state has a total spin of 
S = 10. The most recent DFT parameterization is also 
compatible with INS experiments.^® 


TABLE I. Intra-molecular isotropic exchange interaction pa¬ 
rameters (in meV) as suggested by various authors, compare 
Ref. 16. The spin labels are explained in Fig. 1. 


No. 

Bond (i-j) 

1-6 

1-11 

1-9 

6-9 

7-9 

1-4 

1-3 

1 


(Ref. 16) 

4.6 

1.0 

1.7 

-0.45 

-0.37 

-1.55 

-0.5 

2 


(Ref. 18) 

4.8 

1.37 

1.37 

-0.5 

-0.5 

-1.6 

-0.7 

3 


(Ref. 19) 

5.8 

5.3 

5.3 

0.5 

0.5 

0.7 

0.7 

4 

Jij 

(Ref. 48) 

7.4 

1.72 

1.72 

0 

0 

-1.98 

0 

5 

Jij 

(Ref. 49) 

10.25 10.17 

10.17 

1.98 

1.98 

-0.69 

-0.69 


Figure 2 shows the effective magnetic moment at a 
small external field of H = 0.1 T as a function of tem¬ 
perature. Data of Sessoli®^ and Murrie®® are given by 
symbols. For the theory curves only the Heisenberg part 
of the respective parametrizations is used. Since the 
Heisenberg model is SU(2) symmetric, a FTLM version 
employing total symmetry was used with R = 100 
and Nl = 120 in this case.^® One realizes that the gross 
structure of the magnetic moment (which is proportional 
to VyT) is achieved by all parameterizations especially 
at lower temperatures of T ^ 50 K. A finer inspection 
shows that the maximum is at a too low temperature 
for all parameterizations, so that the experimental low- 
temperature data points are not met. For higher temper¬ 
atures towards room temperature one notices that only 
one parameterization^® (blue curve) closely follows the 
experimental data towards the paramagnetic limit. This 
is not astonishing since this parameterization was fitted 
to the high-temperature tail using a high-temperature 
series expansion. We conjecture that the somewhat too 
large effective magnetic moment of the DFT parame¬ 
terizations, Refs. 16 and 18, at room temperature are 
related to the fact, that these parameterizations con¬ 
tain ferromagnetic interactions whereas a fit to the high- 
temperature behavior, Ref. 19, leads only to antifer¬ 
romagnetic interactions, compare Table 1. In addition 
the antiferromagnetic interactions 1-11 and 1-9 are much 
stronger in Ref. 19. 

Using the recently developed FTLM for anisotropic 
systems^® we could calculate the effective magnetic mo¬ 
ment starting from the DFT parameterization of Ref. 16. 
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FIG. 2. (Color online) Effective magnetic moment of Mni 2 - 
acetate at B = 0.1 T. Data of Sessoli^^ and Murrie®° are given 
by symbols. Observables employing the Heisenberg part of 
parameterizations only are displayed by curves. The param- 
eterizations correspond to those given in Table I from top to 
bottom. 



T (K) 


FIG. 3. (Color online) Effective magnetic moment of Mni 2 - 
acetate at 5 = 0.1 T, same as Fig. 2. The dashed curve shows 
the result of a calculation employing the full Hamiltonian of 
Ref. 16. The powder average is performed over a regular grid 
of 20 directions on the unit sphere.®^ 


Besides the Heisenberg terms of Table I this parameter¬ 
ization contains anisotropic Dzyaloshinskii-Moriya inter¬ 
actions as well as full 3x3 anisotropy tensors for each 
manganese ion. It turns out that the additional terms 
improve the low-temperature data, compare Fig. 3. The 
maximum shifts towards the experimental position and 


the data points for smaller temperatures are much bet¬ 
ter approximated. None of the plain Heisenberg models 
could achieve such an improvement so far. For tempera¬ 
tures above 60 K the anisotropic terms are irrelevant. 


IV. MAGNETIZATION AS FUNCTION OF 
APPLIED FIELD 

The low-temperature magnetization usually provides 
strong fingerprints of the underlying spin-Hamiltonian 
for instance in the case of magnetization steps due to 
ground state level crossings. For Mni 2 -acetate the mag¬ 
netization shows even richer characteristics, since be¬ 
low the blocking temperature a magnetic hysteresis is 
observed.This exciting physical property turns out to 
constitute a problem, when comparing to the theoretical 
equilibrium magnetization. Due to the long relaxation 
times, approximately 2800 hours at T = 2 the ex¬ 
perimental values do not necessarily reflect equilibrium 
values. On the other hand, the theoretical evaluation 
of non-equilibrium observables for a full spin model of 
Mni 2 -acetate is totally out of reach. 


20 

15 
m 
ZL 

2 10 

5 

0 

0 1 2 3 4 5 6 7 

B (T) 

FIG. 4. (Color online) Magnetization of Mni 2 -acetate at T = 
2 K. Color code of curves as above. The dashed curve is 
evaluated for two field values per 1 T field interval only. 

Figure 4 provides two experimental data sets as well 
as various theoretical curves. The data set of Glaser®^ 
was taken on a powder sample whereas the data set 
of Sessoli^^ was taken on a single crystal with a field 
in direction of the tetragonal axis of the S 4 symmetric 
molecule. Both data sets coincide up to H « 2.5 T, 
then the magnetization along the tetragonal axis jumps 
whereas the powder signal smoothly increases with field. 
Already at this point it becomes clear that the measure¬ 
ments cannot reflect equilibrium properties, because the 
magnetization along the tetragonal axis, which is the easy 
axis of this strongly anisotropic molecule,®^ cannot be the 
same as the powder averaged magnetization. 

Interestingly, all theory curves that rest on Heisenberg 
model calculations agree with each other perfectly, which 
is due to the fact that all produce a S' = 10 ground state 
that is largely separated from excited levels. They also 
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agree with the experimental magnetization up to a field 
of B « 1 T. Between 1 T and 2.7 T the experimen¬ 
tal data points stay below the theoretical curves. Above 
B « 2.7 T theory and magnetization along the tetrag¬ 
onal axis meet again. The calculation for the full spin 
model of Mni 2 -acetate as given in Ref. 16 yields a rather 
unexpected result: the powder-averaged magnetization 
stays well below all other theory curves (expected since 
anisotropic), but stays also well below both experimental 
curves (unexpected at least compared to the experimen¬ 
tal powder data). 
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FIG. 5. (Color online) Magnetization of Mni 2 -acetate at T = 
2 K. The dashed curves are evaluated for two field values per 
1 T field interval only. Direction t = (0.0, —0.35682, 0.93417), 
direction p = (—0.35682, 0.93417,0.0). 
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FIG. 6. (Color online) High-held magnetization of Mni 2 - 
acetate at T = 2 K. 


has been measured in megagauss experiments.^^ Interest¬ 
ingly, the magnetization data given in Ref. 54 show pro¬ 
nounced features, likely related to magnetization steps, 
between 180 T and 400 T which could be compatible with 
the parameterization of Ref. 19 (blue curve in Fig. 6). 
As realized already by the authors, this parameteriza¬ 
tion produces a sequence of level crossings between the 
S' = 10 ground manifold and the fully polarized state 
exactly in this field range. 


V. HEAT CAPACITY 


In the following we compare our results with the mea¬ 
surements of Ref. 27 along two different directions, along 
the tetragonal axis and perpendicular to that, i.e. some¬ 
where in the a;j/-plane. Figure 5 presents three theory 
curves: one for the powder average and two along special 
directions. Direction t = (0.0,-0.35682,0.93417) points 
roughly along the tetragonal axis (inclination of about 
20°) and direction p = (—0.35682,0.93417,0.0) lays in 
the xy-p\a,Tie. Both theoretical curves show systemati¬ 
cally larger magnetization values than the experiment. 
This could be for two reasons: either the experimen¬ 
tal curves are not in equilibrium, which is possible at 
T = 2 K where the relaxation time of Mni 2 -acetate is 
long or the parameterization of Ref. 16 is still not yet 
optimal, i.e. the D tensors could be too weak, for in¬ 
stance. Nevertheless, the positive message is, that we 
can now calculate such curves and compare with experi¬ 
mental data. The theoretical costs, by the way, are still 
enormous. For the investigations shown in this article 
about 2 Mio. CPU hours on a supercomputer had to 
be used, since the FTLM procedure has to be completed 
twice for every field value and direction. Therefore, only 
a few field values have been used for the theoretical mag¬ 
netization curves. 

Finally we would like to present the high-held mag¬ 
netization curves. As can be seen in Fig. 6 the vari¬ 
ous parameterizations lead to distinctive differences at 
high helds. The high-held magnetization could be and 


Another observable that was measured very early in 
the history of Mni 2 -acetate is the heat capacity. 

Figure 7 shows the experimental data of Ref. 25 for i? = 0 
(top) and B = 0.3 T (bottom). One notices that the heat 
capacity is rather large and grows steadily with temper¬ 
ature. This is due to a massive contribution from lat¬ 
tice vibrations (phonons) which grows like . There¬ 
fore, heat capacity data of magnetic molecules are usually 
overwhelmed by phonon contributions above T = 5 K. 

This fact becomes obvious when comparing the the¬ 
oretical heat capacity data in Fig. 7 for the Heisenberg 
parameterizations at i? = 0. At low temperatures the 
theoretical values are two or more orders of magnitude 
smaller than the experimental ones. But for the calcula¬ 
tion using the full anisotropic Hamiltonian^® one notices 
that the low-temperature values agree very nicely. We 
think that this is due to a more smeared-out density of 
states at low energies in the anisotropic model whereas 
for Heisenberg systems these levels belong to highly de¬ 
generate multiplets which leads to a different, i.e. much 
smaller heat capacity. 

Interestingly, a magnetic field of H = 0.3 T has a simi¬ 
lar effect. It smears out the density of states due to Zee- 
man splitting. Therefore even for the plain Heisenberg 
models the low-temperature heat capacity increases, but 
still does not agree with the experimental data. For the 
anisotropic spin-Hamiltonian^® the low-temperature heat 
capacity does not change much and still agrees nicely 
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FIG. 7. (Color online) Specific heat of Mni 2 -acetate at B = 0 
(top) and B — 0.3 T (bottom). Data taken from Ref. 25. 
Color code of curves as above. For B = 0.3 T the calculation 
for the anisotropic spin model was averaged over 20 directions. 
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with the experimental data. We conjecture that although 
the energy levels are moved around by the magnetic field, 


the overall structure of the density of states remains very 
similar. 

Summarizing, the specific heat is well reproduced by 
the anisotropic spin-Hamiltonian of Ref. 16 for low- 
temperatures around 1 K. 


VI. SUMMARY AND OUTLOOK 

35 years after its synthesis and 22 years after the first 
measurements^^ of Mni 2 -acetate the Finite-Temperature 
Lanczos Method puts us in a position to evaluate ther¬ 
modynamic functions of really large magnetic molecules. 
It thus complements DFT calculations of such big sys¬ 
tems in so far that one does no longer need to stop half 
way to an understanding of thermodynamic observables. 
The next major necessary step is now to develop tools 
for an evaluation of non-equilibrium properties of such 
big quantum spin systems. 
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